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ABSTRACT 

Several hundred young stars lie in the innermost parsec of our Galaxy. The super-massive black 
hole (SMBH) might capture planets orbiting these stars, and bring them onto nearly radial orbits. 
The same fate might occur to planetary embryos (PEs), i.e. protoplanets born from gravitational 
instabilities in protoplanetary disks. In this paper, we investigate the emission properties of rogue 
planets and PEs in the Galactic center. In particular, we study the effects of photoevaporation, 
caused by the ultraviolet background. Rogue planets can hardly be detected by current or forthcoming 
facilities, unless they are tidally disrupted and accrete onto the SMBH. In contrast, photoevaporation 
of PEs (especially if the PE is being tidally stripped) might lead to a recombination rate as high 
as ~ 10'^® s“^, corresponding to a Brackett-y luminosity L^r-'y ~ 10^^ erg s“^, very similar to the 
observed luminosity of the dusty object G2. We critically discuss the possibility that G2 is a rogue 
PE, and the major uncertainties of this model. 

Subject headings: Galaxy: center - black hole physics - planets and satellites: gaseous planets 


1. INTRODUCTION 

The Galactic center (GC) is one of the most stud¬ 
ied and yet enigmatic places in our Universe. It is ex- 
ceptiona lly crowded: it hosts a supermassive black hole 
(SMBH, iGillessen et al.l 12 009^ . a dense star cluster of 
old stars (iSchodel et al.l 120071 1 . several thous and solar 
masses of molecular, atomic and ionized gas (|Liu et al.l 
I20I2I : iJackson et al.lli99l:lScoville et al.l[2003l l , and a few 


-- - - 

hund r ed young stars (IGenzel et al.ll2003t iPaumard et al.l 

l2?M IBartko et all [2001 iLu et a1.l l2nM 1201 .IH . A rel- 
eyant fraction (^ 20%) of the y oung massiye stars lie 
in the so called c l ockwise disk (IPaumard et al.l I2006I 
IBartko et al.l 120091 : iLu et al.l I2009L 1201311 . a thin disk 
wit h inner radius ^ 0.04 pc and outer radius ^ 0.13 
pc (|Yelda et al.l l2014ll . A small group of B-type stars 
(^ 20 objects, forming the so called S-cluster) orbit 
around the SMBH with semi-major axes < 0.04 pc, with 
isotropically oriented orbital planes, and with high ec- 
cent ricities (approximately following a thermal distribu¬ 
tion, [GiIIiiseVeF2][2QQ93). The formation of the young 
stars in the clockwise disk and in the S-cluster has been 
a puzzle for a long time, because the tidal shear from the 
SMBH disrupts molecular clouds, preyenting star forma¬ 
tion in normal conditions. __ 

A faint dusty object (named G2. IGillessen et aLll2012ll 
has been obseryed on a yery eccentric orbit around the 
SMBH, with a periapse distance o f only ^ 200 AU 
(jWitzel et al.ll20L3 : iPbihl et al.ll2014ll . G2 has transited 
at periapse in Spring 2014, ayoiding complete tidal dis¬ 
ruption. The nature of G2 is still enigmatic: a pure 
gas cloud (IGillessen et al.ll2012l: ISchartmann et al.ll2012l: 


Burkert et al. ~12l: iShcherbakovI 12014 iDe Colle et alJ 

2014 iMcCourt et al.l I2015IL a dust-enshrouded low- 
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mass star (IScoville &: BurkertI 120131 : iBallone et al.l 120131 : 
IWitzel et al.l I2014J1. a low-mass star wi th a proto¬ 
planetary disk (iMurray-Clay fc LoebI 1201211. a star dis¬ 
rupte d by a stellar-mass black hole (iMiralda-Escudd 
I2012I1 . a star t hat underwent partial tidal disruption 
by the SMBH (IGuillochon et all 1201411 . a merger be¬ 
tween two stars iP rodan et al.l 12015), a nd a noya out¬ 
burst (IMeyer fc Meyer-Hofmeistei 20121) haye been pro¬ 


posed to explain its properties Isee lMapelli fc GualandrisI 
120151 for a recent reyiew). 

Planets haye not been detected in the GG so far, 
but the destiny of planets, asteroids and planetesi- 
mals has been inyestigated by seyeral authors. Col¬ 
lisions of planets or asteroids haye been proposed to 
lead to the formation of a dusty torus around SMBHs 
(|Nayakshin et al.l 120121) . The tidal disruption of plan- 
etesimals by the SMBH has been inyoked as mech- 
anism to explain t he daily infrared fl ares of SgrA* 
(iGadez e t al.l 120081: iKostic et alJ 120091 : iZuboyas et al.l 
l2012t iHamers fc Portegies ZwartJ 120151 1. Tidal disrup¬ 
tions of planets are expected to be much less frequent, 
but more dramatic eyents, and might account for the 
possible past actiyi ty of SgrA* ( Reyniytsey et all 12004 
Terrier et al.l 12014 iPonti et alJ 12014 IZuboyas et al l 


20121) . A system composed of a low-mass star and its 


protoplanetary disk is one of the most yiable scenario s 
to explain the G2 object (|Murray-Clay fc LoebI 1201^ . 
Recently, radio-continuum obseryations reyealed 44 par¬ 
tially resolyed compact sources in the innermost ^ 0.1 
pc, interpre ted as candidate photoeva poratiye protoplan¬ 
etary disks (|Yusef-Zadeh et aLll2015D . 

The aim of this paper is to inyestigate the main possi¬ 
ble signatures of planets and planetary embryos (PEs, i.e. 
dense gas clouds produced by lo cal gravitational insta¬ 
bilities in a protoplanetary disk. iKuipeilllO^ iCameroiil 
119781: iBosslllbOTtiDurisen et al.ll2007ll in the GC. In Sec¬ 
tions [2] and [H we describe the mechanisms that might 
produce rogue planets and PEs in the GC, and we esti¬ 
mate the mass loss that planets and PEs undergo because 
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of photoevaporation by the ultraviolet (UV) background. 
In Section m we discuss the observational signatures of 
planets and PEs, with particular attention for the Bry 
line emission), and we suggest that the G2 object might 
be associated with a rogue PE. Section [5] is a summary 
of our main results. 


2. ROGUE PLANETS 


2.1. Tidal capture by the SMBH 

The tidal shear of the SMBH can unbind a planet from 
its star if the initial semi-major axis of the planet orbit 
is 


Op > 19 AU 


0.01 pc 


( m* /4 X 10® M© 

ylOM©/ \ Mbh 


1/3 

( 1 ) 


where d is the periapse of the star orbit around the 
SMBH, m* is the star mass, and Mbh is the SMBH 
mass. Fig. [T] shows Op as a function of d, for m* = 1 
and 10 M©. 

One of the two members of the split binary (generally 
the most massive one) receives a kick that makes it more 
bound to the SMBH, while the other member (generally 
the less massive one) becomes less bound. The less bound 
object might be ejected, while the remaining on e is cap¬ 
tured by the SMBH (e.g. iGinsburg et al.l[^012ll . On the 
other han d, the typic al variation Sy of the velocity of the 
planet is (|Pfah]|l2005ll 



since Sy is lower than the Keplerian velocity around the 
SMBH at d = 0.01 pc (~ 1300 km s“^), it is plausible 
that both the planet and the star remain bound to the 
SMBH. Dynamical simulations are necessary to quantify 
how many planets will be ejected and how many will 
be captured by the SMBH. If the planet is captured, its 
semi-major axis ara.n and eccen tricity Ccap would then be 
pTiiafT^rPerets et al.ll2QM) 


^cap — 1-30 pc 




The tidal split of planets from their stars can be sub¬ 
stantially enhanced by planet-planet scatterings, which 
were shown to be able to scatter Jupite r-mass planets on 
orbit s with semi-major axis> 50 AU (iG h^tterje e et al.l 
120081: iMarois et al.l 120081 iGhatterjee et al.l 1201 Ih . The 
planet could even be eject ed from its i nitial system by 
planet-planet scattering: I Veras et'all (|2009ll estimate 
that ^ 40% of planets in a multiple-planet system are 
ejected by planet-planet scatterings in 2 x 10® Myr. The 
ejection of Jupiter-like planets from their planet systems 
is suppo rted by the obser vation of ‘freely floating’ giant 
planets (iSumi et al.ll2011ll . 

It is even possible that starless planets form directly 
from gra vitational instabilities in accretion disks aroun d 
SMBHs (|Shlosman fc BegelmanlUO^ lNavakshinll200^ . 



Figure 1. Tidal radius for splitting a star-planet system (solid 
black lines) and tidal radius for disruption of the planet (dotted 
red lines) as a function of the star-planet distance from the SMBH. 
Thick and thin solid black lines correspond to a star mass m* = 10 
and 1 Mq, respectively. Thick and thin dotted red lines correspond 
to a mass rrip = 0.02 M© (i.e. a brown dwarf) and 10~® M© (i.e. 
a Jupiter-like planet), respectively. The vertical dashed line is the 
estimated periapse distance of the G2 object. 


Similarly, starless planets could have formed in the same 
star formation episode that gave birth to the clockwise 
disk: according to one of the most popular scenarios, a 
molecular cloud disrupted by the SMBH might have set¬ 
tled into a parsec-scale dense gaseous disk. Gravitational 
instabilities in the gaseous disk might have led to the 
formation of the young massive s t ars that lie in the clock¬ 
wise di sk (e.g. iPaczvnskil 119781 iKolvkh alov fc Snnvaev 
19801 : |Shlosmaii_^Begehnad 19871: IGollin fc Zahn 


1^99; iGammiJ 2001: GoodmanI 120031: PT an fc B lackman 


2005 : iNavakshin et al. 2007: _ Bonnell fc Rice 20081 


MaDelli e t al.ll200 8r Hobbs fc Navak shiDll200 9HAlig et all 

20111 IMa, cellj et all 120121 120131 lEucas et all 1201310 


and also to the formation of starless giant planets 
(|Shlosman fc Begelmanl[l989D . In such case, the initial 
orbit of the planets would be inside the clockwise disk, 
and then gravitational interactions with stars or other 
planets might plunge the planets on a mor e radial orbit. 
For example, iMurrav-Glav fc LoebI (|2012| ) predict that 
^ 1/100 of the entire population of low-mass objects in 
the clockwise disk could have been delivered onto highly 
eccentric orbits by two-body interactions with massive 
stars. 

What happens to a planet that is captured by the 
SMBH on a very eccentric orbit? Tidal disruption is 
unlikely, as the planet should have a radius Tp larger 
than 


rt = 1.3xl0"%m- -^ — 

0.01 pc V Mbh ) VIO" 


(5) 

i.e. ^ 2000 (d/0.01 pc) times the radius of Jupiter. Fig.[T] 
shows the behavior of rt as a function of d, for a planet 
mass rrip = 10“® M© and for a brown dwarf of mass 
rUp = 2 X 10“^ M©. 
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2 . 2 . Photoevaporation 

The planet will suffer continuous atmosphere evapora¬ 
tion by the UV field of the voung stars i n the central 
parsec. According to iMurrav-Clav et aP (1200911 (using 
the simplified eq. 19), the mass loss rate by atmosphere 
evaporation is 

. eTTrlLuy/iAirD'^) 

ivjMC — /^ / ’ 

Grup/rp 

where rp is the planet radius, Luv 10 "^° erg s" 
the ionizing luminosity of massive stars in the GC, 
the fraction of the UV luminosity that goes into heat 
(~ 4'l4>o — 1 ) where (j) is the average energy of ionizing 
photons, and (j)Q ~ 13.6 eV is the ionization energy for 
atomic H), and D is the distance of the massive stars 
from the rogue planet. We assume that D ~ 0.1 pc, 
since this is the outer rim of the clockwise disk in the 
GG (lYelda et al.ll 20 Tl . 

For large values of Cp, eq. ([ 6 ]) implies that the num¬ 
ber of ionizations is larger than the number of ionizing 
photons reaching the planet surface. 


^ _ 7rr2Luv/(4 7rZl2) 

(/<o/(l - e) ■ 


(7) 


rate can be estimated as 


R = 




-3 


A 2 3 1 ^max . 2 3 a f ^-l- ™prot \ 

= 4 TT ttB n_^_ Tp In- ~ 4 TT aB Tp In I - - — 1 

fp \ Ph J 


dr 


2/3 



= 4 TT ttB n~p r® 

(6) 

- 3 X 10®®s"® 

^ is 


e is 



In 


(ri+ mprot) 
Ph 


9.7 


( 10 ) 


where aB — 2-6 x 10 ^^cm^s ^ is the case B re¬ 
combination coefficient for Hydrogen (at a temperature 
^ 10^ K), and we have chosen Cmax as the radius where 
the density of the evaporated gas drops to the den¬ 
sity of the hot medium: ruprot ?^(?'max) = Ph (the nor¬ 
malization of the logarithmic term is appropriate for 
n+ = 10^cm“^, Ph = 10“^^gcm“^). The corresponding 
emission measure (EM) is EM= f n^dV ~ 10"^® cm“^, 
for Tp = 10 ^° cm and = 10 ^ cm“^. 

Equation cni) assumes that the evaporated gas is 
nearly completely ionized. This is a good approxima¬ 
tion, because the recombination time scale 


Therefore, we employ an evaporation rate 


M = min(MMC, Wprot Qp) = 


^p Guv UT-prot 

4772 ^ 


()>o/mprot \ 

- 1^,(1 - e) 


(Gmp)/rp 




^uv 


\ lO^^ergs" 


D 


-2 


X mm 


0.98 e 


(to5Iu)"( 


10 ^° cm 


O.lpc 


,( 8 ) 


where mprot — 1-67 x 10 24 g is the proton mass. The 
bottom panel of Fig. [2] (dotted line) shows the behavior 
of M as a function of the planet radius. 

The average kinetic energy of the evaporating ions is of 
the order of ~ <p~4>o = c ['i^o/(l ~ e)]; implying a velocity 
of the order of Ug ~ 30kms“^ (for e ~ 0.3). This velocity 
is generally of the same order of magnitude as the escape 
velocity Vesc = v^( 2 GTOp)/rp. 

Thus, the density of the ionized gas that evaporates 
from the star is 


ivi 7 —Q / jri \ 

n_L = -— = 10 cm - 77 ^-r 

4 TTmprot r2'Cg y6xl04°gs 

^30kms“^^ ^ 10 ^°cm^^ 

The central panel of Fig. [5] (dotted line) shows n_|_ (ob¬ 
tained assuming Vg = Vesc = \/ ( 2 Gmp)/rp), as a func¬ 
tion of the planet radius. 

If the velocity of the evaporating matter is close to 
the escape velocity from the planet, the density profile 
becomes n(r) ^ n+(r/rp)“^7^, and the recombination 


Gec(r) = -^-ru- ~ 3.8xl0'^s ( , 3 ) 

anrij-ir] \ 10^ cm / 


aBn+(r) 


is longer than the ionization time scale 



1 

cr(^)Luv/(47rZl2^) 


~ 1800s 


/104°ergs-i\ / D 
\ Guv / V 0.1 pc 
( 12 ) 


where we assumed that the the ionization cross_section of 
neutral Hydrogen is cr{cj)) ~ 6.3 x 10“^® cm^{(j)/(j)o)~^ ~ 
2.1 X 10“^®cm2, jg appropriate for e = 0.3, i.e. 

<^ = </,o/0.7~ 19.4 eV. 

Furthermore, in this case, we can ignore the possibility 
of shocks with the hot medium, because the relatively low 
values of Vg often lead to a Mach number M = Vg/cs ^ 1, 
where Cs is the sound speed (slow-wind case). In the Ap¬ 
pendix we discuss the case in which A4 >> 1 (fast- 
wind case). Photoevaporation in the fast-wind approx¬ 
imation leads to a recombination rate about one order 
of magnitude lower with respect to the slow-wind ap¬ 
proximation, but shocks occurring in the fast-wind case 
enhance the recombination rate by a factor depending on 
M (see Appendix rA)l . 

In the following, we refer to CASE 1 (see Table 1) as 
the model where M, n_|_ and R have been calculated from 
equations E m and [ini respectively (i.e. R has been 
calculated in the slow-wind approximation, with Vg = 

y '2 Gmp/rp and Guv = 10“^° erg s“^). Fig. [ 2 ] (dotted 
line) shows M, and R in CASE 1, as a function of 
the planet radius. 


3. ROGUE PROTOPLANETARY EMBRYOS 
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Table 1 

Summary of model properties. 


Name 

M (gs-i) 

n_(. (cm ^) 

R (s-i 

) 

Tidal stripping 

CASE 1 

from eq. 

from eq. 

from eq. 


no 

CASE 2 

from eg. Il5l 

from eg. Il4| 

from eq. 

10 

no 

CASE 3 

- 

- 

from eq. 

[2^ 

yes 


Note. — Columns 2, 3 and 4 specify how M, n_|_ and R were calculated in each 
model. CASE 1 corresponds to a photoevaporating planet (eq.[^ Section 2) in the slow- 
wind approximation (eg, llUl , with rrip = 10~^ ^0> '^g ~ ^^2 G mp/rp and r^uV ~ 
10^^ erg s~^. CASE 2 corresponds to a photoevaporating cloud (eq. Section 3) 
in the slow-wind approximation (eq. with mp = 10~^ = ca = 10 km 

s~^ and r^uV ~ 10“^^ erg s~^. CASE 3 corresponds to a photoevaporating cloud 
undergoing tidal stripping (rp > r^), with vg = cs = 10 km s~^ and r<uv = 10^^ 
erg s~^. If the cloud is being stripped, we do not need to derive M and n_j. in order 
to calculate R, since R = (see eq. ESI and the discussion in the text). 



Figure 2. From top to bottom: recombination rate associated 
with photoevaporation (R), evaporating gas density (n+), and 
mass loss rate by photoevaporation (M) as a function of the planet 
(or PE) radius. Black dotted line: CASE 1 (photoevaporating 
planet, Table 1). CAS E 1 is valid onl y for rp < 2 x 10^^ cm (i.e. 
for T < 1, see Section l3.ll and eg. 1131 for details). Red solid line: 
CASE 2 (photoevaporating cloud, Table 1). The two blue dashed 
lines in the top panel are the minimum and maximum value of the 
recombination rate measured for the G2 cloud in the Br 7 line since 

One of the two main competing scenario^ for the for¬ 
mation of Jupiter-like planets and brown dwarfs predicts 

^ In this paper, we neglect the competing core accretion model 
(i.e. the accretion of a gaseous atmosphere on a rocky core , 
e.g. IWetherilllllQSg : [Bodenheimer &; Pollackin-9861: ILissauenflOCTl , 

which will be considered in forthcoming studies. 


that they form through local gravitation al instabilities 
in the outer parts of a protoplanetary d isk (iKuipeiillQSlI : 
iBossUlQQ^IlQQSal lbl nHelled et al.ll2008D . The local insta¬ 
bility might produce a protoplanetary embryo (PE): a 
coreless gas clump with density ~ 10 “® g cm“^ (much 
lower than typical planetary and stellar densities), and 
large radius (corresponding to a Jeans length Aj ~ 
few AU). Then, these gas clumps cool down, contract¬ 
ing to radii and densities typical of brown dwarfs or gi- 
ant gaseous planets. R ecent population synthesis studies 
(|Forgan fc Ricell2013ll indicate that PEs have a mass of 
a few to tens of Jupiter masses, and that a large fraction 
of PEs (~ 40 — 90 %) end up forming brown dwarfs (with 
mass 0.02 Mq). 

One of the main predictions of this theory is that PEs 
can form only in the outer parts of proto planetary disks 
(^ 1 0 — 50 AU distance fr om the star , IPorean fc Ricel 
[Ml. or even > 100 AU, IM)^: IRafikovI I200l 7 

where the self-gravity of gas is sufficiently strong with 
respect to stellar gravity and radiation pressure. Thus, 
a PE is an excellent candidate for being tidally captured 
by the SMBH, before it can contract to a Jupiter-like 
configuration. Furthermore, starless PEs might form di¬ 
rectly from gravitational instabilities in a dense gaseous 
disk surrounding the SMBH (see previous section). 


3.1. Photoevaporation 

For a radius of the PE < rt -- 1.3 X 

10^^ cm (d/0.01 pc) (see eq. [5]), the PE will avoid tidal 
disruption during the capture by the SMBH. From eq. |H1 
we infer that the mass loss by evaporation for a PE 
is M ^ 3.5 X 10^® g s“^ for rp = 5 x 10^^ cm (and 
e = 0.3), corresponding to n+ ~ 2.9 x 10® cm“® for 
Vg = ^2Gmplrp. Thus, the recombination rate for an 
evaporating PE would he R ^ 3.1 x 10^^ s“^ (using eq. 
dni), corresponding to an EM ^ 10®^ cm ®. However, 
this result neglects the optical depth encountered by the 
ionizing photons before reaching the PE surface. 


r = 


5 X 10 


-4 


n+ 


.107 


cm 


cr(0) n+ 

r)^(- 


-3/2 


trecis ) 


dr ~ 


/lO^o ergs-7 


10^° cm/ \ Luv 


D 


0.1 pc 


where we approximated the neutral fraction in the evap¬ 
orating gas at radius r as tion/Uec(?') (which is correct as 
long as tion < Uec(?•))• 

If we use the CASE 1 approximations at 
rp = 5 X 10^^ cm, we get r ~ 200, implying 

that the estimate of R is too high, since the 
ionizing photons would be unable to photoevap- 
orate the surface of the PE. In general, our ap¬ 
proximations break down when r 1 , i.e. when 
rp lAU(n+/107cm“®)“^ (Luv/ 10 ‘*°ergs“ 7 ) (T)/0.Ipc)“ 
in particular, CASE 1 holds only if rp 2 x lO^^ cm. 

Therefore, it is more appropriate to derive the mass 
loss rate M with the assumption that the PE is a pure 
gas cloud (such as a proto planetary disk). Following 
IMurrav-Clav fc LoebI (|2012ll . the mass loss rate of a gas 
cloud with radius rp, which is undergoing photoevapora¬ 
tion, can be expressed as M ^ 4 tt r^ mprot''^■- 1 -Cs, where 
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the number density of the photoevaporated ions (n+) can 
be calculated assuming balance between recombinations 
and photoionizations at the base of the wind, 


between the bo w shock of the stellar wind and the 
hot m edium is (|Burkert et al.l 120121 : iScoville fc Burkerll 

[Mil 


n+Ri (Luv/(4 7r0i:)^))^/2 (aer-p) (14) 


where (j) is the average energy of ionizing photons. Thus, 
the mass loss rate is 





^prot 


Csr 


3/2 

P 




- 6.7 X lO^^gs"^ 


) ( 

lOkms V 


rp 

lO^^cm/ 


Tuv 


10 "^° ergs" 



and the optical depth due to the evaporating gas becomes 
r ^ 0.5, independent of all the considered parameters, 
apart from 0. This result shows that eq. (HI can be 
applied for a PE, since the absorption in the photoevap- 
orative wind does not stop the photoevaporation itself. 

Furthermore, if we substitute eq. da into eq. da, 
the recombination rate becomes 




47rr^Luv , 

- - —^ In 

4:7rD‘^(j) 


^max 




(16) 


where Qp is the number of ionizing photons reaching the 
PE surface in the optically thin approximation (eq. [7]). 
We have R > Qp because r-max ^ Cp, so that the number 
of available ionizing photons is much larger than Qp. 

The central and bottom panel of Fig. [2] (solid line) show 
the behavior of n+ and M, as derived from equations 1141 
andfTSl respectively. Finally, the top panel of Fig. [7] (solid 
line) shows the recombination rate R obtained combining 
eq. [14] with eq. [10] (or, equivalently, using eq. [T6t . In the 
following, we define this model as CASE 2 (see Table 1). 


3.2. Tidal stripping enhancement 
The tidal radius rt of a PE can become similar (or 
smaller) than its radius rp for d ^ 0.01 pc (Fig. [T]); in 
such case, mass loss by tidal stripping might become non- 
negligible at some point in the orbit. When rt < rp, 
we can estimate the mass-lo ss rate by tidal stripping as 
(jMurrav-Clav fc Loebll20l3l 


Mtid 



( "l 

f A7bh \ 

V lO^kms-1 J 

\Ax lO6M0y 


-1/3 

(17) 


where we assume that the surface density of the PE is 
p(rp) ~ TOp/(dirrp) (appropriate in the case of an isother¬ 
mal density profile), a nd v \ is the radial compon ent of 
the orbital velocity Up ([Murrav-Clav fc Loebll201^ . 

Which is the fate of the stripped material? How could 
we observe it? The tidally stripped material might un¬ 
dergo shocks with the high-temperature medium. The 
stagnation radius Vg where ram pressure is balanced 


M Vg 
4 TT Ph lip 


1/2 


= 2 X 10^® cm 


M, 


tid 


8.7 X 1022gs-i 


f W f 

Vl 0 kms“i/ \ 


10-21 gcm-3y/^ /1300 km s"^ 
Ph y V Up 


1/2 


(18) 


where Vp 1300kms i) is the Kepler velocity of a 
planet orbiting the SMBH at 0.01 pc, and ph = lO”^! 
g cm”® is the density of the h ot medium from X-ray 
measurements ([Yuan et al.ll2003[l . 

Thus, the maximum luminosity that can be emitted by 
the tidally stripped material in such shocks is 


Tmax — TT 


\ ^prot / 


1 


1 


1 (-) ^prot ^p — Q -^tid '^g '^p 


,36^^^ ^-1 


~ 1.4 X 10 ergs 

( 


Mfid 


8.7 X 1022 gs- 


VlOkms 0 (1300kms i)^ 


where we used equations m and m for rg and Mtid, 
respectively. This luminosity is very high but can be 
achieved only if: (i) the stripped material reaches the 
stagnation radius rg (but we estimate that the stripped 
gas needs t ^ 300 yr to reach the stagnation radius, rg fv 
IQi^cm, if it travels at Vg ^ lOkms-i); (ii) all of the 
kinetic energy is converted into (observable) luminosity. 
Thus, we expect that the luminosity due to these shocks 
is much lower than Tmax, and we will neglect it in the 
rest of this paper. 

On the other hand, the stripped material will be 
expos ed to photoevaporation b y the massive young 
stars (|Murrav-Clav &: Loebll201^ . The presence of the 
stripped material modifies the geometry of the PE and 
m ight strongly enhance phot oevaporation (as discussed 
in iMurrav-Clav fc Loebll2?)T^ . To correctly evaluate the 
new geometry is beyond the aims of the current paper, 
but we can account for this correction in the following 
way. The surface of the stripped material will take part 
in photoevaporation, and will produce an approximately 
spherical wind. 

The recombination rate in the wind will be of the same 
order of magnitude as the rate of ionizing photons that 
can be absorbed by the tidally stripped material, that is 


^ 47rr2^, Luv 
47ri42())g 


( 20 ) 


where rgtr is the maximum radius reached by the tidally 
stripped material. We evaluate rstr as 

^str ^ T ft 'Cg, (^4) 

i.e. as the radius of the PE, plus the distance travelled by 
the stripped material (moving at velocity Vg) during the 
time tt between the instant when the PE radius becomes 
equal to the tidal radius (rt = rp), and the time of the 
periapse passage. 
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Figure 3. Bottom panel: radius rgtr (eq. 1211 assuming periapse 
distance ~ 200 AU, eccentricity e = 0.976 and velocity Vg = 10 
km s“^), as a function of PE radius. A value r^tr = Tp means 
that Cp is smaller than rt for the entire orbit (in this case, CASE 3 
is the same as CASE 2). Top panel: recombination rate R in 
CASE 3 (Ta ble 1), as a function of PE radius. R has been derived 
from eq. 1201 assuming that R = Qtid- In both panels, magenta 
solid line: PE mass rrip = 0.02 Mq; green dot-dashed line: PE 
mass nip = 10“® Mq. The two blue dashed lines in the top panel 
are the minimum and maximum value of the recomb ination rate 
measured for the G2 cloud in the Bry line since 2004 JPfuhl et al.l 

EoH). 

Fig. [3] shows the recombination rat43 R (assumed to 
be equal to Qtid) and the radius Tstr, as derived from 
equations [20l and [2T1 respectively. In Fig. [3l we assume 
Vg = 10 km s“^, periapse distance = 200 AU and eccen¬ 
tricity e = 0.976, i.e. the same periapse and eccentricity 
as the G2 object. If rp > rt along the entire orbit, we 
assume R = O.STprb (where T^rh is the orbital period). 
In principle, this kind of calculation can be applied also 
to planets, but for the orbit of the G2 cloud rt is always 
larger than rp if rp iS 1.5 x 10^^ cm (for rrip = 10“^ Mq). 
In the following, we refer to the model presented in Fig. [3] 
as GASE 3 (see Table 1). 

4. DISCUSSION 

4.1. Luminosity of rogue planets in the GC 

Do we have any chances of detecting rogue planets or 
PEs in the GG with current or forthcoming facilities? 
The K (2.1 p,m) and L' (3.8 fj,m) magnitudes of a rogue 
planet at the distance of the GG (~ 8 kpc) are tok 32 

The calculations of IMurrav-Clav fc LoebI II2012I 1 are somewhat 
simil ar to ours, even if they apply to a protoplanetary disk. How¬ 
ever, IMurrav-Clav &: Locbf (120121) do no t ba lance ionizations and 
recombinations. As illustrated by eq. eg, such balance is not 
required in the case of planets or PEs, because rmax ^ How¬ 
ever, when the tidal stripping enhancement is important, the ra¬ 
tio Tmax/’^str should be of the order of unity, and R Qtid- 
For example, in the case of disks with radius ^ 10 AU producing 
ph otoevaporative winds with n ~ 10 ^cm“^, the equations used 
by IMurrav-Clav Sz LoeO (120121) imply a recombination rate much 
larger than the ionization rate, and the material would become 
neutral in ~ 5 days. 



Figure 4. Br 7 luminosity (Z/Bry) as a function of the planet (or 
PE) radius. The lines and colors are the same as used in Figs. [2] 
and [3] In particular, dotted black line; CASE 1; solid red line: 
CASE 2; green dash-dotted line: CASE 3 with rrip = 10“^ M 0 . 
The two blue dashed lines are the minimum and maximum val ue 
of Z/Bry for the G2 cloud observed since 2004 (jPfuhl et al.ll2Ql^ . 

and Tot/ 26 for a mass rup ^ 10“^ Mq (lAllard et al.l 
1200112001 . respectively. Therefore, such rogue planet 
would be invisible for current and forthcoming facili¬ 
ties (30-m class telescopes are expected to observe stars 
down to Tok 24). On the other hand, some processes 
might take place that enhance the chances of observing a 
planet, such as tidal disruption, atmosphere evaporation 
and bow shocks. In the previous section, we analyzed 
these processes for both planets and PEs. 

Fig. m shows the luminosity of the Bry line (2.166 
/rm) derived from Lbt'j = 2.35 x 10“^^ ergs“^ R/ae, 
where the recombination rate R was calculated in the 
previous section. The most optimistic prediction for 
a Jupiter-like planet (mp = 10“^ Mq and rp = 10^° 
cm) is a Bry luminosity ~ 5 x 10^® erg s“^, if 

the planet atmosphere is undergoing photoevaporation. 
At the distance of the GC, this corresponds to a flux 
10 “^^ergs“^ cm“^, which is far below the sensitivity 
of existing and forthcoming instruments (a 24-hr integra¬ 
tion with VLT SINFONI can theoretically detect a Bry 
flux ^ 5 X 10“^® ergs“^ cm“^ with a S/N of ^ 10). 

Thus, the only chance of detecting a rogue planet in the 
GC is that it is disrupted and accretes onto the SMBH. 
If a portion of the planet mass is accreted by the SMBH, 
this might lead to a flare with bolometric luminosit y 
< 2 X erg s“^ (from eq. 45 of iZubovas et al.ll2012f) . 
This event is very unlikely, because the tidal disruption 
occurs only if the periapse distance from the SMBH is 
1.6 AU (rp/10^°cm). We expect a planet disruption 
rate Rdis ^ 10 ~^/p yr“^ (where /p is the number of 
planets per star in the GC, and 10~® vr~^ is the tida l 
disruption rate of stars estimated bv IAlexandeill 200 '^ . 
This process is very rare, but might be relevant for ex¬ 
plaining flares in other galactic nuclei. 

4.2. Luminosity of rogue PEs in the GC 
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Theoretical models (e.g. iWuchte^ 119991: iBossI 120051 : 
iHelled et al.ll2006t iHelled fc Bodenheimerll2010ll suggest 
that a 0.001 Mq PE does not collapse immediately (be¬ 
cause it is optically thick) and needs to cool for 1 Myr 
at a bolometric luminosity ~ 10“® Lq ~ 10^^ergs“^. At 
the distance of the GC, this luminosity corresponds to 
tol' .^25, which is far below the limits of current and 
forthcoming observational facilities. 

The photoevaporation of the ‘atmosphere’ of a PE 
(with a radius of rp = 5 x 10^^ cm) could lead to a Bry 
luminosity of ^ 5 x lO^^erg s“^, considering CASE 2 es¬ 
timates for the mass loss by photoevaporation, assuming 
Luv = 10“^° erg s“^. If the PE is undergoing tidal strip¬ 
ping, this luminosity might be boosted by more than 2 
orders of magnitude. A LBr 7 ~ 10^° erg s“^ can be ob¬ 
served with current 8m telescopes, and is not far from 
the actual luminosity of the G2 cloud (^ 6 x 10^° erg 



Tidal disruption of a PE occurs if the distance of the 
PE from the SMBH is d ~ 750 AU (rp/5 x lO^^cm), sug¬ 
gesting that the cross section for PE disruption is a factor 
of ~ 10® larger than the cross section for planet disrup¬ 
tion (neglecting gravitational focusing). This leads to a 
disruption rate of ^ /pe yr~^, where /pe is the number 
of PEs per each star. 

How many PEs can form in the GC and for how long do 
they survive? This question contains a number of ma¬ 
jor uncertainties and we can just suggest a few hints. 
PEs are quite elusive objects and are expected to be 
relatively short-lived before they contract to Jupiter- 
like size (^ iWuchterl 119991 : IHelled et al.ll200d 

iForgan fc Ricell2?)T^ . but in the GC they could be effi¬ 
ciently separated from their parent star, and the strong 
UV flux might substantially slow down their cooling (and 
collapse). Furthermore, we could reverse our question, 
and use the non-detection of TBr 7 ~ 10®^ erg objects 
to constrain the frequency of PEs in the GC. 

4.3. Rogue PEs and the G2 cloud 

The G2 cloud is the only observed object (so far) that 
shares similar properties with a photoevaporating and 
partially stripped PE. G2 probably originated from the 
clockwise disk in the G C, and has a very h igh eccentric¬ 
ity, e ^ 0.976 ± 0.007 (iPfuhl et al.ll2014ll . Both these 
constraints on the orbit are fairly consistent with the hy¬ 
pothesis of tidal capture of a PE (initially formed in the 
clockwise disk) by the SMBH, or of a dynamical ejection 
of the PE from the clockwise disk. It is now clear that 
G2 survived its peri apse passage at a d istance of ^ 200 
AU from the SMBH (iWitzel et al.ll201^ . but some of its 
material was tidally stripped. A PE with radius ^ 10^^ 
cm is tidally stripped if the periapse is ^ 160 AU: such 
PE would not be completely tidally disrupted at a dis¬ 
tance of ^ 200 AU from the SMBH, but it would suffer 
some tidal stripping. 

Fig. 0] shows that the Brq luminosity of a PE (as 
derived from photoevaporation and partial tidal strip¬ 
ping) matches the one observed for G2, if the PE has 
nip Ri 10“® Mq and rp r; 2 x 10^® cm. If the main trigger 
of the Bry emission is photoevaporation by the intense 
UV background in the GC, we expect the Bry luminosity 
to remain roughly constant during the PE orbit. This is 
fairly consistent with the fact that the Bry luminosity 



Figure 5. Uncertainties on the Br -7 luminosity of a PE. Red 
filled cross-hatched area: CASE 2 by varying Ljjy from 10^® erg 
s“^ (lower bound) to erg s~^ (upper bound). Dot-dashed 

black line: same as CASE 2 with Ljjy = 10'^® erg s~^ but in 
the fas t-wind approximation (with Vg = 10 Cg, see Appendix m 
eg.lAlll. Solid black line: Br -7 luminosity from shocks in the fast - 
wind approximation (with Vg = 10 Cg, see Appendix [A] eq. IA2|I . 
Green filled area: CASE 3 with Luv = 10“^^ erg s~^, by varying 
the PE mass from 10“^ Mq (upper bound) to ~ 0.2 Mq (lower 
bound). Dashed green lines: rup = 10“^ Mq (top), 0.02 Mq 
(intermediate), 0.1 Mq (bottom). Horizon tal blue filled are a: the 
observed Br -7 luminosity of the G2 cloud IjPfuhl et al.ll2Qi^ b 

of G2 remained nearly constant over 10 yr t|Pfuhl et al.l 

120111 . 

G2 was detect ed also in the L ' band (ttil' ^ 14, 
iPfuhl et 'all 1201411 . iGillessen et al.l (|2012h suggest that 
the continuum L' emission comes from small (~ 20 nm), 
transiently heated dust grains with a total warm dust 
mass of ~ 2 X 10^® g. The grains mig ht be warmed up by 
an iiiner source (e.g. a low-mass star. iScoville fc BurkertI 

120131 : [Ballone et al.ll2013l : [Witzel et al.ll2014li 7 or by some 

external mechanism (UV heating, shocks, etc.). In the 
PE scenario, the minim um PE radius neces sary to reach 
Ll' 2.1 X 10®®ergs“^ (IWitzel et al.ll201‘3 l is 


5.5 X lO^^cm 


Al' 


2.1 X 10®®ergs“ 


1/2 


Adust 
560 K 


m) 


where Adust 560 K is the estimated du st temperature 
(from L' — M' ^ 0.3. IGillessen et ^l2012f) . We note that 
■Cmin strongly depend on Aiust and that the estimate of 
Adust is very uncertaiifl. 

While the scenario of a wind-enshrouded low-mass 
star would natural ly explain the continuum L' emission 
(|Witzel et al.ll20fll . we cannot reject the hypothesis that 
a PE, embedded in the hot dense medium of the GC, 
might host sufficient warm dust to power the observed 
L' emission. Thus, a rogue PE might be a vi able scenario 
to ex plain G2 and other G2-like objects (e.g. IPfuhl et'ffi] 
120141 suggest that the object G1 is related to G2). 


4.4. Discussion of uncertainties 


^ It should be noted that neither 
the radiation absorbed by the PE ('^ 
10 ^^ ergs“^ [Luv/( 10 ^'^ ergs“^)] [D/( 0.1 pc)]“^[rp/( 10 ^^ cm)]^), 
nor the energy released in shocks around it (< 2.5 x 

10^0 ergs“^ [ph/( 10 “^^ gcm“^)] [rp/( 10 ^^ cm)]^ [up/( 1000 kms“^)]^) 
can balance with Ll/, unless rp ^ 5 x 10^^ cm. 
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The analytic model discussed in this paper relies on 
the possibility that planets form in the GC. The GC 
might be a hostile environment for planet formation, 
given the high UV background, the high temperature 
and the strong tidal field by the SMBH. Discussing the 
chances that planets form in the GG is beyond the aims of 
this paper. Here, we just mention that planetesimals and 
asteroids were re cently investigated as possible sources 
of SMBH flares (iCadez et al.ll2008t iKostic et al. 120091 : 
iZubovas et al.l l2012t iHamers fc Portegies ZwartI 120151 1 . 
and there are even some observational hints for the ex- 
iste nce of protoplanetary d isks in the innermost par¬ 
sec (lYusef-Zadeh et ^l2015f) . Further observational evi¬ 
dence is necessary, to confirm that planets and planetary 
objects form in the GG. 

Planets cannot be directly observed with current and 
forthcoming facilities. Only photoevaporating PEs are 
sufficiently bright to be detected. PEs are theoretically 
predicted objects but have not yet been observed. They 
can form only if gravitational instabilities in a proto¬ 
planetary disk lead to fragmentation, and to the forma¬ 
tion of self-gravitating clumps. Recent work has shown 
that gravitational instabilities can lead to the formation 
of PEs only in the o utermost reg ions of a protoplane¬ 
tary disk (> 100 AU, iBolevll^OQr) . where self-gravity is 
stronger. This enhances the probability that PEs become 
unbound with respect to their initial system, but might 
imply that their formation is endangered by the strong 
UV field in the GG. Furthermore, even if hydrodynam- 
ical simulations show that clumps can form, they can¬ 
not predict self-consistently whether these clumps sur¬ 
vive further evolution. 

Recent estimates suggest that PEs remain in the 
pre-collapse phase (the one considered in this paper) 
for ^ 10^ yr and for ^ 10'^ yr i f their mass is ^ 
10~^ and 10~^ Mr?^, respec tively (|Helled et al.l 120061 : 
iHelled fc Bodenheimerl 120101 ). Erom an observational 
perspective, while the detection of massive giant plan¬ 
ets ( ^ 10“^ Mq) at distance >>10 AU from the cen- 
tral star might favor the gravitational instability scenari o 
(jMarois et al.l 120081 : see Figure 3 of iPeoe et al.l [20T^ . 
there is no strong evidence supporting the existence of 
PEs. Since PEs are such elusive objects, it is quite diffi¬ 
cult to quantify the uncertainties in our model. 

Besides the debate on the very existence of PEs, Eig. [5] 
shows the impact of the main sources of uncertainties on 
the predicted Br-y luminosity of a photoevaporating PE. 
In particular, we focus on (1) the adopted approximation 
for the wind (fast or slow-wind approximation), (2) the 
flux of the UV background (from 10^® erg s“^ to 10"^® erg 
s“^ in the innermost 0.1 pc), (3) the mass of the PE. 

Eig. [5] shows that the Br-y luminosity of a photoevap¬ 
orating PE is uncertain by several orders of magnitude. 
We stress that the recombination rate (and thus the Br- 
y luminosity) does not depend on the mass of the PE 
in CASE 2. The mass of the PE is important only if 
we assume that the PE can be tidally stripped by the 
SMBH (CASE 3), because rt depends on the PE mass. In 
the fast-wind approximation (i.e. the gas is ejected with 
Vg >> Vesc, as discussed in detail in Appendix lAT) , the 
Br-y luminosity due to photoevaporation (from eq. [AH) 
is lower by about one order of magnitude with respect to 
the slow-wind approximation, but the gas may undergo 


shocks with the hot medium, and this enhances the Br-y 
luminosity to a value (from eq. IA2I) comparable with (or 
even higher than) the slow-wind approximation. 

5. SUMMARY 

In this paper, we investigated the possible observa¬ 
tional signatures of planets and planetary embryos (PEs) 
in the GC. If planets and PEs form in the central parsec, 
several mechanisms can separate them from their par¬ 
ent star (e.g. tidal shear by the SMBH or planet-planet 
scatterings), and bring them onto a very eccentric orbit 
around the SMBH. It is even possible that starless PEs 
and planets form directly from gravitational instabilities 
in a dense gas disk around the SMBH (such as the one 
that might have given birth to the clockwise disk of young 
massive stars). 

We have shown that both planets and PEs suffer from 
photoevaporation (Eig. [5]) due to the intense UV back¬ 
ground. The emission measure associated with such pro¬ 
cess is relatively low for planets (EM^ 10^® cm“®) and 
much higher, although very uncertain, for PEs (EM« 
j^q 50-56 gg^gg of PEs, tidal stripping can en¬ 

hance the effect of photoevaporation, leading to an even 
higher EM (up to EMr; 10®° cm“®, Fig. |3|). This means 
that a photoevaporating PE with radius ~ 5 x 10^® cm 
might reach a Br-y luminosity LBr 7 ~ 5 x 10^® erg s“^ if 
it is not tidally stripped, and TBr 7 ~ 5 x 10®^ erg s“^ if 
it is partially tidally stripped (Eig. SJ. This value, while 
uncertain, is remarkably similar to the observed Br-y lu¬ 
minosity of the G2 dusty object. Furthermore, a PE with 
radius ^ 5 x 10®® cm can emit the same Ll' luminosity as 
G2, if it contains the sufficient amount of dust at temper¬ 
ature ~ 600 K. In our model, the L' luminosity is emitted 
from a smaller area than the Br-y line, since the former 
is due to dust inside the PE, while the latter is produced 
by the photoevaporative wind. This can account for the 
fact that the observed L' emission is mor e compact than 
the emission in Br-y (jWitzel et al.ll20fll . If G2 is a PE 
with radius ~ 5 x 10®® cm, we expect its lifetime to be 
tiife ~ 10® yr (top/10“®Mq) (5 x 10®'’g s“®/M), but tidal 
stripping can reduce tufe significantly (down to ~ 100 

yr)- 

Our results are affected by several uncertainties. First, 
PEs are theoretically predicted objects, but elusive to 
observe. Their properties and survival time are uncer¬ 
tain. Furthermore, the luminosity of photoevaporating 
PEs strongly depends on several quantities (e.g. PE 
mass, UV background luminosity, wind speed), as shown 
in Fig. O In a follow-up work we will investigate the hy- 
drodynamical evolution of PE models embedded in a UV 
background. Furthermore, the frequency of PEs in the 
GG, and the probability that they are captured by the 
SMBH deserves further study. Our preliminary results 
open a new exciting window on GC’s environment. 
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APPENDIX 

A. THE FAST-WIND APPROXIMATION 

In the main text, we discussed the case in which 
Vg ^ Vesc- While this is the more likely scenario, in this 
appendix we also consider the case in which the initial 
velocity Vg >> Vesc- If Tg >> Vesc, we can assume that 
Vg remains approximately constant, so that the gas den¬ 
sity scales as n(r) ^ n+ (r/rp)“^. Thus, we can estimate 
the recombination rate as as 

dr 

= 47raB A - 

where as ^ 2.6 x 10“^^ cm^ s“^ is the Case B radia¬ 
tive recombination coefficient for H (at a temperature 
of ^ 10^ K), and Cmax (typically ^ Cp) is the outer 
limit of the wind-dominated density profile (either the 
radius where the interaction between the wind and the 
high-temperature medium produces a shock - see be¬ 
low -, or the radius where the wind density drops below 
that of the surrounding medium). A recombination rate 
R ^ 1.2 X 10^^s“^ corresponds to an emission measure 
EM=JnldV - 3 X lO^® cm-^. 

The wind that evaporates from the planet will also 
undergo a shock with the high-temperature medium in 
the GC. Equation [18] in Section [32] provides the stagna¬ 
tion radius rg where ram pressure is balanced between 
the bow shock of the stellar wind and the ho t medium 
([Burkert et al.ll2M^ iScoville &: Burker'tll2013ll . 

We estimate the combined effect of the shock and the 
p hotoevaporatio n with a simplified version of the results 
of iDvsonf (jl975[ l: we assume that the results along the 
direction of motion are approximately valid for all direc¬ 
tions, obtaining a recombination rate 


i? ~ 4 TT aB Tp 




2.5 X lO-'^s 


, 33-1 / '^+ 


)\- 


Vl0^cm“3/ VlO^^cm 


(A2) 


Equation IA2I holds only if A4^ >> 1, where M = Vg/cg 
is the Mach number, a nd C s is the sound speed (the nor¬ 
malization used in eq. IA2I adopts A4 = 5, corresponding 
to Vg = 50 km s“^, and Cg = 10 km s“^). Thus, the 
contribution of shocks increases the recombination rate 
by a factor of ~ 8 for A4 = 5. In Fig. [SJ we compare the 
Br-y luminosity of a PE in the slow-wind approximation 
and in the fast-wind approximation. 
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